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The effect of Mg codoping and Mg deposition mode on the Fe distribution in (Ga,Fe)N layers 
grown by metalorganic vapor phase epitaxy is investigated. Both homogeneously- and digitally-Mg 
codoped samples are considered and contrasted to the case of (Ga,Fe)N layers obtained without 
any codoping by shallow impurities. The structural analysis of the layers by high-resolution 
transmission electron microscopy and by high-resolution- and synchrotron x-ray diffraction gives 
evidence of the fact that in the case of homogenous- Mg doping, Mg and Fe competitively occupy 
the Ga-substitutional cation sites, reducing the efficiency of Fe incorporation. Accordingly, the 
character of the magnetization is modified from ferromagnetic-like in the non-codoped films to 
paramagnetic in the case of homogeneous Mg codoping. The findings are discussed vis-a-vis 
theoretical results obtained by ab initio computations, showing only a weak effect of codoping 
on the pairing energy of two Fe cations in bulk GaN. However, according to these computations, 
codoping reverses the sign of the paring energy of Fe cations at the Ga-rich surface, substantiating 
the view that the Fe aggregation occurs at the growth surface, fn contrast to the homogenous 
deposition mode, the digital one is found to remarkably promote the aggregation of the magnetic 
ions. The Fe-rich nanocrystals formed in this way are distributed non-uniformly, giving reason for 
the observed deviation from a standard superparamagnetic behavior. 

PACS numbers: 68.55.Nq, 75.75.Cd, 81.15.Gh,75.50 Pp, 61.10 Nz 



I. INTRODUCTION 



One of the principal goals of nowadays materials sci- 
ence is to develop tools allowing for the bottom-up assem- 
bling of nanostructures with predefined properties and 
functionalities. It was theoretically suggested^ that the 
codoping of magnetic semiconductors with shallow impu- 
rities affects the self-assembly of magnetic nanocrystals 
during epitaxy, and therefore modifies both the global 
and local magnetic behavior of the material. This con- 
cept was qualitatively corroborated by experimental data 
for (Zn,Cr)Te (Ref.0) and (Ga,Fe)N (Refs.ii). In the 
latter case it was proven that by doping (Ga,Fe)N with 
donors (Si) the charge state of the Fe ions is reduced 
from Fe 3+ to Fe 2+ , affecting the binding energy of the 
Fe pairs, and suppressing the ions' aggregation even at a 
high Fe precursor flow rate^ The quenching of Fe aggre- 
gation, and therefore of either crystallographic or chemi- 
cal phase-separation, results in an overall reduced magne- 
tization of the layers, particularly at high temperatures, 
where the paramagnetic response of dilute spins is small. 

While the effect of doping of (Ga,Fe)N with donors 
may be considered experimentally clarified, the role of 
intentional doping with acceptors (Mg) rises several re- 
markable issues, namely: (i) does the introduction of ac- 
ceptors in a feasible concentration modify the Fe charge 
state to Fe 4+ or do the holes remain bound to the Mg 



ions?— £ (ii) in either case, does the introduction of accep- 
tors affect, similarly to donors, 3 ' 4 the aggregation of Fe? 
(hi) which is the effect of Mg-codoping on the magnetiza- 
tion of the layers? (iv) is the process of aggregation essen- 
tially a bulk effect [as in annealed (Ga,Mn)As] (Refs.0, 
8) or it occurs at the growth surface, as suggested by 
our preliminary experimental 3 - and theoretical 9 results 
for (Ga,Fe)N? (v) can the chemistry and spatial distri- 
bution of Fe-rich nanocrystals be controlled by codoping 
protocols? 

In the present work, we report on the effect of Mg dop- 
ing and deposition mode, considering three different se- 
ries of (Ga,Fe)N samples: (i) without any codoping; (ii) 
with Mg continuously provided during the growth pro- 
cess (homogeneous codoping) and (hi) with digital Mg- 
codoping ((5Mg)ji2. Our x-ray diffraction (XRD), trans- 
mission electron microscopy (TEM), and superconduct- 
ing quantum interference device (SQUID) data show that 
homogeneous Mg codoping reduces the Fe incorporation 
with the consequent suppression of Fe-rich nanocrystals 
formation, as preliminary reported^ Interestingly, ut- 
terly different effects are found in the case of digital 
codoping. In this case, the self-assembling of Fe-rich 
nanocrystals is promoted, and their chemical composi- 
tion and spatial distribution affected. Thus, the present 
analysis complements, by providing fundamental infor- 
mation on the role of the mode of acceptors incorpo- 
ration into the system, the extensive work other o 5 ! 11 ^ 



and w tf 15 ' 16 have carried out on controlling the struc- 
tural and magnetic properties of (Ga,Fe)N. The data pre- 
sented here demonstrate that various modes of codoping 
may serve, in a bottom-up fashion, to hamper or to pro- 
mote phase separation and the self-assembling of Fe-rich 
nanocrystals in GaN. 

These experimental results are confronted with our 
ab initio calculations carried out according to a method- 
ology developed and employed earlier— to confirm the 
experimental observation 3 that the aggregation of Fe in 
GaN occurs at the growth surface. Other relevant first 
principles studies include works on the electronic states 
brought about by Ga-substitutional Fe ions in GaN4*il 
It was found that the Fe 3+ ion, in addition to acting as 
a mid-gap acceptor, gives rise to a hole trap level local- 
ized at 0.3 eV above the top of the valence band,— sub- 
stantiating earlier theoretical suggestions^ and experi- 
mental findings^ Recently, the energetic and magnetic 
properties of Fe incorporated on a clean and Ga-bilayer 
GaN(OOOl) surface have been evaluated^ within a for- 
malism similar to oursi* It was found that a Fe atom oc- 
cupies preferentially, by ^2.0 eV, a surface substitutional 
Ga site rather than a near-surface interstitial site. More- 
over, the authors conclude that under N-rich conditions 
the Fe atoms incorporate into Ga substitutional sites at 
the top bilayer with no tendency to migrate towards the 
bulk. 



Here, we present computed values of cohesive energies 
for various Fe, Mg, and also Si clusters in the bulk and 
at the Ga-terminated GaN(OOOl) surface, assuming that 
the Fe ions occupy Ga-substitutional positions, in agree- 
ment with the first principles studies referred to above»i2 
The computational data allow us to interpret the effect 
of codoping by either Mg or Si on the incorporation of 
the Fe ions. We also put forward an explanation for the 
highly non-random distribution of Fe-rich nanocrystals 
in the film volume, and relate it to the magnetic proper- 
ties of this ferromagnetic/semiconductor nanocomposite 
system. 

The manuscript is structured as follows: the experi- 
mental details on the fabrication and on the characteri- 
zation methods are given in Sec. II, where a table listing 
the samples considered and their significant parameters is 
included. In Sec. Ill we discuss the structural and mag- 
netic analysis of the two differently Mg-doped series of 
(Ga,Fe)N epilayers. We then present, in Sec. IV, the re- 
sults of first principles computations carried out within a 
model developed to account for the Fe-Mg interplay ob- 
served experimentally and to compare it with the Fe-Si 
case. A series of often unanticipated conclusions emerg- 
ing from the comparison of experimental and computa- 
tional results are outlined in Sec. V. Within the final sec- 
tion, we provide a summary and an outlook prompted by 
the results. 



II. SAMPLES AND EXPERIMENTAL 
PROCEDURE 

A. Samples 

Reference (Ga,Fe)N and two series of (Ga,Fe)N lay- 
ers codoped with Mg have been fabricated by metalor- 
ganic vapor phase epitaxy (MOVPE) in a 200RF AIX- 
TRON horizontal system at a substrate temperature 
T S = 850°C with trimethylgallium (TMGa), ammonia 
(NH 3 ), ferrocene (Cp2Fe) and biscyclopentadienyl mag- 
nesium (Cp2Mg) as precursors. The 1/im thick doped 
layers are deposited onto a I /im thick GaN buffer grown 
on c-plane sapphire (AI2O3) substrates following the pro- 
cedure described elsewhere^ employing a TMGa pre- 
cursor flow rate of 12 standard cubic centimeters per 
minute (seem). For the homogeneously Mg-doped lay- 
ers TMGa and both doping sources Cp2Mg and Cp2Fe 
are opened simultaneously during the layer growth. The 
(Ga,Fe)N:<5Mg layers are produced by growing 60 peri- 
ods of (Ga,Fe)N layers between Mg thin layers*^ The 
enrichment of the sample surface with nitrogen is en- 
sured during the entire growth procedure by keeping the 
NH3 flow rate constant. The Cp2Mg source flow rate is 
350 seem, giving - according to secondary ion mass spec- 
troscopy (SIMS) measurements - a mean Mg content of 
«10 19 cm~ 3 (iM g = 0.02%) in the <5-doped samples and 
a concentration ten times higher in the homogenously 
doped ones. The Fe concentration is varied by changing 
the Cp2Fe source flow rate from 50 to 450 seem. A sum- 
mary of the main parameters employed to fabricate the 
samples considered in this work is presented in TablelU 



B. Secondary-ion mass spectroscopy 

The total Fe and Mg content in the layers has been 
determined via secondary-ion mass spectroscopy (SIMS). 
An undoped GaN sample implanted with 10 16 cm~ 2 of 
56 Fe at 270 keV is used as reference for calibration of the 
system. 15 In the (Ga,Fe)N reference layers the total Fe 
concentration in the films is found to increase with the 
nominal Cp2Fe flow rate. 



C. X-ray diffraction 

Structural analysis of the samples is performed by 
high- resolution x-ray diffraction (HRXRD), and syn- 
chrotron x-ray diffraction (SXRD). Conventional x-ray 
diffraction measurements are carried out on a X'Pert Pro 
material research diffractometer at an energy of 8 keV 
with the following configuration: a hybrid monochroma- 
tor and a 0.25° divergence slit are placed in front of the 
incoming beam and a PixCel detector with 1 mm active 
length is used, in order to reduce the diffuse scattering 
around the substrate peak and to improve the resolution. 



TABLE I: (Ga,Fe)N references and Mg-doped (Ga,Fe)N sam- 
ples with the corresponding Fe and Mg mean concentrations 
as evaluated via SIMS. The doping method is indicated with 
h for Mg-homogeneous, and 5 for Mg-digital doping. 



sample 


type 


Cp 2 Fe 


Cp 2 Mg 


XFc 


^Mg 






[seem] 


[seem] 


[%] 


[%] 


S907 


ref. 


100 





0.06 





S911 


ref. 


200 





0.11 





S910 


ref. 


300 





0.21 





S994 


S 


100 


350 


0.15 


0.02 


S995 


S 


200 


350 


0.20 


0.016 


S993 


S 


300 


350 


0.38 


0.016 


S1193 


S 


50 


350 


0.08 


0.03 


S1192 


ref. 


50 





0.08 





S1406 


ref. 


300 





0.20 





S1407 


h 





350 


0.00 


0.15 


S1402 


h 


50 


350 


0.03 


0.23 


S1405 


h 


100 


350 


0.04 


0.22 


S1404 


h 


200 


350 


0.07 


0.20 


S1403 


h 


300 


350 


0.11 


0.19 


S1415 


h 


400 


350 


0.12 


0.17 


S1416 


h 


450 


350 


0.15 


0.17 


S1427 


h 


300 


250 


0.12 


0.07 


S1428 


h 


300 


150 


0.10 


0.01 



Due to the high intensity beam and to the enhanced res- 
olution obtained with the mentioned configuration the 
diffraction peaks originating from embedded nanocrys- 
tals are resolvable. However, to obtain an enhanced sig- 
nal to noise ratio, SXRD measurements have been per- 
formed at the BM20 Rossendorf Beamline (ROBL) of 
the European Synchrotron Research Facility in Greno- 
ble, France. The samples are measured in coplanar ge- 
ometry using a photon energy of 10.005 keV. Diffraction 
scans are acquired along the growth direction [001], and 
in a range between the (002) diffractions of GaN and the 



(006) of A1 2 3 ^ 



D. Transmission electron microscopy 

The high-resolution transmission electron microscopy 
(HRTEM) studies are performed on cross-sectional sam- 
ples prepared by standard mechanical polishing followed 
by Ar + ion milling at 4 kV for about 1 h. Conventional 
diffraction contrast images in bright-field imaging mode 
and high-resolution phase contrast pictures have been ob- 
tained from a JEOL 2011 Fast TEM microscope operat- 
ing at 200 kV and capable of an ultimate point-to-point 
resolution of 0.19 nm and allowing to image lattice fringes 
with a 0.14nm resolution. Additionally, energy disper- 
sive spectroscopy (EDS) analysis via an Oxford Inca EDS 
equipped with a silicon detector provides information on 



the local composition. Selected area electron diffraction 
(SAED) and fast Fourier transform (FFT) procedures 
are employed to study the scattering orders and the d- 
spacing for respectively the larger and the smaller em- 
bedded nanocrystals. 



E. SQUID magnetometry 



The magnetic moment of the samples has been mea- 
sured in a Quantum Design MPMS XL 5 SQUID magne- 
tometer between 1.85 and 400 K and up to 5 T following 
the methodology described elsewhere . 20 ' 21 

As demonstrated previously^! 6 - ' 15 ' 16 the magnetic re- 
sponse of (Ga,Fe)N layers can be considered as repre- 
sentative for semiconductors containing transition metal 
ions at a concentration near the solubility limit. In par- 
ticular, it contains a sizable diamagnetic component from 
the sapphire substrate, which is carefully subtracted ac- 
cording to the procedure detailed recently^ 20 ' 21 The re- 
maining magnetic moment m(H) is: (i) purely param- 
agnetic and assigned to substitutional Fe 3+ ions if the 
system is dilute, while (ii) it is build up from the param- 
agnetic part as well as from a ferromagnetic-like compo- 
nent and from an antiferromagnetic contribution linear 
in H - these latter, both stemming from various mag- 
netically ordered nanocrystals with high Fe content^ 6 - - 
if the material undergoes phase-separation. 

Due to the fact that the Fe-rich nanocrystals are 
characterized by relatively high ordering temperatures, 
a significant variation of m(H) at low temperatures 
allows us to assess quite precisely the concentration 
of paramagnetic ions a; para . Following the proce- 
dure employed previously, 6,1 - we obtain x para by fit- 
ting gfiBSxparaNoABsiH, T) to the difference between 
the experimental values of the magnetization measured 
at 1.85 and 5.0 K, where AB S (H,T) is the difference 
of the corresponding paramagnetic Brillouin functions 
ABs(H,T) = B s (H,lMK)-B s (H,5K). We adopt here 
parameters suitable for Fe 3+ ions in GaN, i.e., the 
spin 5=5/2, the corresponding Lande factor g = 2.0, 
and the cation density iVo =4.4 x 10 22 cm~ 3 , treating 
a^para as the only fitting parameter. Its value is then 
employed to calculate the paramagnetic contribution 
at any other temperature according to M(T, H) = 
g/iBSxp^^NoBs^H^T), which is then subtracted from 
the experimental data to obtain the magnitude of the 
magnetization M{ erTO (H, T) and its saturation value from 
which we establish a;f crro , a lower limit of the Fe concen- 
tration in the Fe-rich nanocrystals. We assume the mag- 
netic moment per Fe ion to be Sg/j,B =2.0/xb, the value 
determined for e-FesNJ^ Since, in general, the distribu- 
tion of ferromagnetic nanocrystals is highly non-uniform, 
a^fcrro is intended as a mean value. 



F. Electron spin resonance 



Electron spin resonance (ESR) spectra have been ob- 
tained with a Bruker Elexsys E 500 spectrometer op- 
erating in the microwave X-band (~ 9.8 GHz) equipped 
with a pumped helium cryostat. The measurements have 
been performed at a temperature of 2 K at the microwave 
power of f=a 2 mW. The magnetic held is applied parallel 
to the crystal c-axis. To allow for lock-in detection, its 
magnitude is modulated at 100 kHz with an amplitude of 
5 Oe. The integrated amplitude of the Fe 3+ ESR signal, 
omitting the region in the immediate vicinity of residual 
resonances corresponding to g = 2.0, serves as a relative 
measure of the concentration of Fe 3+ ions J£ 
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III. EXPERIMENTAL RESULTS 

We discuss here the results of the SIMS, XRD, TEM, 
SQUID, and ESR measurements carried out according 
to the methodologies specified in the previous section. 
The findings provide information on the Fe incorpora- 
tion, distribution, concentration, and on the correspond- 
ing magnetic properties of the considered samples. In 
order to demonstrate that we can control the Fe in- 
corporation and distribution by selecting an appropri- 
ate codoping mode, the experimental results for three 
series of samples, namely: (Ga,Fe)N (reference sam- 
ples), (Ga,Fe)N:Mg (homogeneous Mg-codoping), and 
(Ga,Fe)N:<5Mg (digital Mg-codoping) are compared. 



A. Fe incorporation — SIMS data 

In Fig.Q] the Fe concentration as obtained from SIMS 
data averaged over the layer thickness is reported. As 
seen, the two employed Mg codoping modes have a con- 
siderable, but reverse effect on the Fe incorporation, 
that is on the total Fe concentration in the layers for 
a given Fe precursor flow rate. In particular, homoge- 
neous Mg-codoping reduces substantially the Fe concen- 
tration, which is - on the other hand - significantly en- 
hanced in case of digital Mg-codoping in comparison to 
the (Ga,Fe)N reference films. It is interesting to estab- 
lish how the effects described above evolve with the Mg 
concentration. According to the SIMS data presented in 
Fig. [2] only above a certain threshold of the Mg concen- 
tration, namely for iMg > 0.03%, Mg appears to prevail 
in thwarting the efficient incorporation of the magnetic 
ions. 



B. Fe-rich nanocrystals - XRD and TEM data 

As it can be seen in Fig. [3] in (Ga,Fe)N films, grown 
at a sufficiently high Fe precursor flow rate of 300 seem, 
there are signatures of crystallographic phase separa- 
tion witnessed by the HRXRD peaks corresponding to 



FIG. 1: Total Fe content as evaluated from SIMS measure- 
ments as a function of the Fe-source flow-rate for the reference 
(Ga,Fe)N and for the correspondent (Ga,Fe)N:Mg x Mg = 0.2% 
layers, as determined from SIMS. 
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FIG. 2: Mg concentration as determined by SIMS vs. Cp2Fe 
source flow for (Ga,Fe)N:Mg (squares) and (Ga,Fe)N:5Mg 
(triangles). GaN:Mg reference samples [homogeneous (<) and 
digital (t>)] are included. All layers with iMg < 0.03% show 
Fe crystallographic phase separation for Cp2Fe source flows 
as low as 50 seem. 



the expected e-Fe3N secondary phase, detected also by 
HRTEM, as displayed in the inset to Fig.[3l In contrast, 
in the correspondent (Ga,Fe)N:Mg sample grown with 
the same parameter, but homogeneously doped with Mg, 
the Fe ions are less efficiently incorporated and the sam- 
ple is dilute, as confirmed in Fig.[3]by the absence of both 
additional diffraction peaks in the HRXRD spectra and 
of nanocrystals in the HRTEM images. 

In Fig.|H the HRXRD spectrum for a reference dilute 
(Ga,Fe)N sample grown at the relatively low Fe precur- 
sor flow rate of 100 seem, is compared to those obtained 
for a homogeneously- and for a digitally-Mg codoped 
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FIG. 3: HRXRD curves for a reference (Ga,Fe)N grown at 
Fe precursor flow rate of 300 seem with £-Fe3N embedded 
nanocrystals and for the corresponding dilute (Ga,Fe)N:Mg 
layer with £ M g = 0.2%. Insets: HRTEM images - (a): e-Fe 3 N 
embedded in the (Ga,Fe)N reference sample; (b) no precipi- 
tates in the corresponding (Ga,Fe)N:Mg layer. 



layer grown under the same Fe flow rate. While for 
the (Ga,Fe)N reference and for the homogeneously Mg 
codoped (Ga,Fe)N:Mg layer there is no evidence of crys- 
tallographic phase separation and only the diffraction 
peaks from the GaN buffer and from the sapphire sub- 
strate are detected, the (Ga,Fe)N:<5Mg sample gives clear 
evidence of diffraction arising from secondary phases. 
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FIG. 4: HRXRD spectra for a reference dilute (Ga,Fe)N sam- 
ple and for corresponding layers homogeneously and digitally 
codoped by Mg, respectively. All films were grown at the Fe 
precursor flux rate of 100 seem. 

In order to identify the detected Fe-rich phases, SXRD 
measurements have been carried out on a series of 
(Ga,Fe)N:<5Mg samples fabricated with constant nominal 
<5Mg-doping and by varying the nominal Fe concentra- 



tion throughout the series. In Fig. [5] the SXRD spectra 
for these samples are reported together with the corre- 
sponding (Ga,Fe)N phase-separated (nominal Fe content 
above 0.2%) reference. Several diffraction peaks close 
to the AI2O3 (006) peak are observed and are labeled 
from 1 to 8 in the inset to Fig.[SJ The Fe-rich phases in 
the (Ga,Fe)N phase-separated sample are all observed in 
the Mg-codoped ones, however their density is higher in 
the <5Mg-doped layers as inferred from the higher peak 
intensity and the narrower full-width-at-half-maximum 
(FWHM) of the peaks from the nanocrystals in the <5Mg- 
doping case compared with those in (Ga,Fe)N. The iden- 
tified Fe-rich phases, their d-spacing and their size, as 
obtained directly from the FWHM of the SXRD peaks^ 
are summarized in Table[IIJ Some of the phases identi- 
fied, like e.g. the galfenol Fei-^Ga^ ones, are not present 
in phase-separated (Ga,Fe)N when not <5Mg-doped. 
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FIG. 5: SXRD spectra of highly Fe-doped (Ga,Fe)N (*) 
and (Ga,Fe)N:5Mg samples with different Fe concentrations 
(0.15%<x<0.38%) giving evidence of the presence of sec- 
ondary Fe-rich phases in the layers. Inset: detail of diffraction 
peaks from the secondary phases - gaussian fit of the peaks. 
The phases corresponding to the diffraction peaks labelled 1 
to 8 are listed in TablelTTl 

The FFT of the Moire fringes in the HRTEM images 
of selected nanocrystals, reveals the presence of mainly 
e-Fe3N in the (Ga,Fe)N reference and a variety of Fe- 
rich phases with d-spacings as reported from SXRD in 
the (Ga,Fe)N:JMg ones. The EDS spectra show that the 
Fe signal is significantly enhanced in the region of the 
nanocrystals, confirming their Fe-rich nature. 

Another important issue concerns the spatial distribu- 
tion of the Fe-rich nanocrystals. SIMS measurements on 
the <5Mg-doped samples evidence an inhomogeneous dis- 
tribution of Fe in the layers along the growth direction, 
as seen in Fig.lSJa). Up to 300 nm below the surface 
the Fe concentration reaches its maximum for all <5Mg- 
doped samples. Additionally, in the sample containing 
0.38% of Fe, an increased concentration of the magnetic 



TABLE II: Diffraction peaks from Fe-rich phases in 
(Ga,Fe)N:JMg, as observed in Fig.0 



peak 


d-spacing 


av. size 


phase (hkl) 


Nr. 


(nm) 


(nm) 




1 


0.2188 


13 


£-Fe 3 N (002) 


2 


0.2102 


43 


£-Fe 2 N (Oil) 


3 


0.2077 


60 


e-Fe 3 N (111) 


4 


0.2057 


64 


Fe 3 Ga (220) 


5 


0.2048 


38 


£-Fe 2 . 4 N (221) 


6 


0.2020 


25 


a-Fe (110) 


7 


0.1998 


25 


Feo. 7 Ga .3 (101) 


8 


0.1914 


24 


7'-Fe 4 N (200) 



ions is resolved close to the interface. This is also ob- 
served in TEM images, where in the highly <5Mg-doped 
layer the nanocrystals are found 300 nm below the sam- 
ple surface - as seen in Fig.[6jb) - in two narrow re- 
gions separated by 180 nm and close to the interface - 
as evidenced in Fig.[5Jc). In the undoped (Ga,Fe)N sam- 
ple, the nanocrystals lie around 200 nm below the sample 
surface.—^ 




0.5 1.0 
Depth (pm) 

FIG. 6: (a) SIMS profiles of (Ga,Fe)N:<5Mg proving the in- 
homogeneous distribution of Fe in the layers. (b),(c): TEM 
bright-field images of nanocrystals in the GaN:Fe:<5Mg sample 
S993 (xf c = 0.38%) - (b) close to the sample surface; (c) at 
the interface GaN-buffer/(Ga,Fe)N:<5Mg-layer. 



C. Magnetic properties and Fe concentrations — 
SQUID data 

A strict correlation between the presence of Fe- 
rich nanocrystals and the emergence of high tempera- 
ture ferromagnetism has already been demonstrated for 
(Ga,Fe)N (Ref .iGJGJ) . In fact, the SQUID data can 
serve to quantify the concentrations of Fe ions contribut- 



ing to ferromagnetic and paramagnetic signals, respec- 
tively, as described in Sec. lHEl 

Since the concentrations of Fe and Mg are similar, the 
question arises about a contribution of the Mg impuri- 
ties to the magnetic response. From ESR measurements, 
we can state that no paramagnetic signal originates from 
holes bound to Mg acceptors in GaN:Mg, presumably be- 
cause of compensation by hydrogen centers. Likely due to 
the same reason, the paramagnetic signal is - within 3% - 
independent of the orientation of the c-axis of the samples 
with respect to the magnetic field, indicating that the Fe 
ions retain the 3+ charge state also in (Ga,Fe)N:Mg. 
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FIG. 7: (Color online) Mean concentrations of Fe ions 
contributing to ferromagnetic response as a function of Fe 
flux rate, as measured by SQUID for (Ga,Fe)N (circles); 
(Ga,Fe)N:Mg (squares); (Ga,Fe)N:5Mg (triangles). Dashed 
lines are guides for the eye. 

In Figs.[7]and[8]the amount of Fe contributing to the 
ferromagnetic and paramagnetic signals, is respectively 
given. We see that the onset of ferromagnetism, which 
occurs at about 100 seem of the Fe-precursor-flow-rate in 
the case of (Ga,Fe)N, is shifted to up to 300 seem and 
down to 50 seem for homogeneous and digital Mg codop- 
ing, respectively. Thus, the present data supports the 
assignment of ferromagnetism to the presence of Fe-rich 
nanocrystals revealed by XRD and TEM, as discussed 
in the previous subsection. Moreover, they demonstrate 
directly how a particular mode of codoping can serve to 
control high temperature ferromagnetism in (Ga,Fe)N. 

As seen by comparing Figs. [7] and [3 in the region where 
ferromagnetism sets in, the concentration of paramag- 
netic Fe ions saturates and then decreases. It is also 
clear that homogeneous Mg-codoping reduces not only 
the ferromagnetic contribution but also the magnitude of 
the paramagnetic signal, in accordance with the induced 
quenching of the efficiency of Fe-incorporation. Interest- 
ingly, the data obtained earlier for (Ga,Fe)N:Si (Ref.y), 
and displayed also in Fig.El show that Si codoping, in 
contrast to Mg codoping , enhances the concentration of 
paramagnetic Fe ions, at least at low Fe flow rates. 
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FIG. 8: (Color online) Mean concentrations of Fe ions con- 
tributing to the paramagnetic response as a function of the 
Fe flow rate, as measured by SQUID for (Ga,Fe)N (circles); 
(Ga,Fe)N:Mg (squares); (Ga,Fe)N:5Mg (triangles). Results 
obtained earlier for (Ga,Fe)N:Si (Ref.y) are shown for com- 
parison (diamonds). Dashed lines are guides for the eye. 



In Fig.[9j the total Fe concentration from SIMS is com- 
pared to the sum of ferromagnetic and paramagnetic 
components. It is seen that the concentration from SIMS 
is systematically higher in the ferromagnetic region, in- 
dicating that some Fe ions are not captured by the em- 
ployed analysis of the SQUID data in this regime. We 
assign the missing contribution to the presence of antifer- 
romagnetic Fe-rich nanocrystals stabilized by the hostJ^ 



The clear correlation between the appearance of 
ferromagnetic-like features and the presence of Fe-rich 
nanocrystals referred to above substantiates the fact that 
Fe aggregation accounts for the high temperature ferro- 
magnetism of (Ga,Fe)N. As shown in Fig.flOl the fer- 
romagnetic component of the magnetization Mf(H) of 
(Ga,Fe)N exhibits the features characteristic for high 
temperature ferromagnetic semiconductors. In particu- 
lar, for any orientation of the magnetic field, the satura- 
tion magnetization is much larger than the spontaneous 
one. However, the spontaneous magnetization persists 
to above room temperature, indicating that for some 
nanocrystals the blocking temperature is higher than 
300 K. Actually, the shape of the normalized M(H), par- 
ticularly the presence of a small but non-zero sponta- 
neous magnetization even at 300 K, indicates that for 
some nanocrystals the superparamagnetic limit has not 
been reached, whereas for others M{ is described by a 
temperature independent Langcvin-typc function: 23 



M f /M f Sat = tacnh(H/H ), 



(1) 



where Hq is a sample dependent parameter, evaluated 
here to be 800 ± 200 Oe. 

It has been suggested 2 - 3 that dipole-dipole interactions 
between densely packed magnetic constituents may lead 
to a magnetization that can be accurately parameter- 
ized by a temperature independent Langevin-type func- 
tion (Eq.[l]), with AttHq > M f Sat constituting the mag- 
nitude of the local magnetization. In particular, these 
interactions may result in a domain or vortex structure 
within individual magnetic nanocrystals, that reduces 
their magnetization in weak magnetic fields. Further- 
more, as evidenced by the TEM and SIMS measurements 
discussed in the previous subsection and illustrated in 
Fig. E] the distribution of magnetic nanocrystals is highly 
non-uniform, promoting their dipole-dipole interactions. 
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FIG. 9: (Color online) Total concentrations of Fe ions deter- 
mined by SIMS (full symbols) and obtained as a sum of the 
ferromagnetic and paramagnetic contributions determined by 
SQUID magnetometry (open symbols) for (Ga,Fe)N (circles); 
(Ga,Fe)N:Mg (squares); (Ga,Fe)N:c5Mg (triangles). Dashed 
and dotted lines are guides for the eye. 



D. Concentration of Fe 3+ ions — ESR data 

In order to asses the microscopic origin of the para- 
magnetic component we have carried out ESR measure- 
ments. As in the previous studies^ we observe a spec- 
trum that is assigned to non-interacting Fe 3+ ions in a 
wurtzite crystal field. As shown in Fig.[l2l the Fe concen- 
tration determined from the ESR signal intensity scales 
consistently in both (Ga,Fe)N and (Ga,Fe)N:<5Mg with 
the concentration of paramagnetic ions drawn from the 
SQUID data. At the same time, however, the ESR signal 
is dramatically reduced in (Ga,Fe)N:Mg. 

The integrated peak amplitude of the Fe 3+ magnetic 
transitions at the magnetic field parallel to the crystal c- 
axis serves as a relative measure for the concentration of 
these isolated paramagnetic ions, 15 calculated according 
to the following procedure: since the spectra are collected 
by means of a lock-in technique, they are proportional to 
the first derivative of the absorption with respect to the 
magnetic field. 
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FIG. 10: (Color online) Normalized magnetization at vari- 
ous temperatures of sample S995 [(Ga,Fe)N:<5Mg, xf c = 0.2%] 
with embedded nanocrystals. Points - experimental data; 
solid line - fit a Langevin-type function with Ho — 800 Oe. 
Arrow: sweep direction of the magnetic field. 



First a linear background is subtracted from the spec- 
tra which is afterwards numerically integrated to yield a 
signal proportional to the absorption. These integrated 
spectra usually show a broad background, which is re- 
moved by substracting a polynomial of an order between 
4 and 9. Then the peaks of interest are integrated again 
to get their integrated intensity, the integration ranges 
being chosen to be the same for all samples. These inte- 
grated peak areas are proportional to the total number of 
spins and are therefore normalized by the sample volume. 

To convert these integrated peak areas into absolute 
concentration values, they are multiplied by a calibration 
factor (the same for all samples), which is chosen such 
that it minimizes the sum over all samples of the squares 
of the differences between the Fe 3+ concentration values 
obtained by ESR and SQUID according to the procedure 
described in Sec. Ill El The formula that is minimized is 
E = J2 (alESRi — #SQUiD.) j where a is the calibration 

i 

factor, Ii is the integrated peak intensity for sample i 
obtained from ESR and Xi is the Fe concentration in 
sample i as obtained by SQUID. The minimization leads 



to a 



J2 ^ESR^SQUID, / I X)^E 



'ESR, 



The Fe transition at g m 2 has not been employed for 
quantification, since it overlaps with a signal stemming 
from impurities in the substrate and varying in intensity 
from wafer to wafer. 



rates our statement, that most of the present Fe is indeed 
isolated and in the 3+ valence state. By comparing the 
ESR spectra of (Ga,Fe)N with those of (Ga,Fe)N:Mg, it 
is observed that the concentration of Fe 3+ is generally 
decreased when homogeneously codoping (Ga,Fe)N with 
Mr 
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FIG. 11: (Color online) ESR for (Ga,Fe)N:Mg (squares) and 
reference (Ga,Fe)N (circles). Upper panel - one of the se- 
lected Fe 3+ transition and double integral for the samples 
grown with 200 seem Cp2Fe flow rate: (Ga,Fe)N (circles, 
solid line) and a (Ga,Fe)N:Mg (squares, dashed line), respec- 
tively. Lower panel - Fe 3+ concentration, as determined by 
ESR (hlled symbols) and SQUID (open symbols). 



2. (Ga,Fe)N:SMg 

Also for the digitally Mg-doped layers a significant con- 
sistency between the SQUID and ESR values for Fe 3+ is 
found and summarized in Fig. 1121 In the <5Mg doped sam- 
ples, the concentrations of Fe 3+ are systematically and 
notably lower than the overall Fe concentration as deter- 
mined by SIMS, in agreement with the observation that a 
substantial amount of Fe is incorporated in secondary Fe- 
rich phases. Furthermore, we generally observe that the 
Fe 3+ concentration as determined by ESR is lower in the 
digitally doped samples than in the reference (Ga,Fe)N. 



1. (Ga,Fe)N:Mg 

As seen in Fig. llll the concentration of paramagnetic 
isolated Fe 3+ as determined by ESR is consistent with 
the values obtained by fitting the SQUID magnetisation 
curves according to Sec. Ill El This agreement corrabo- 



IV. COMPUTATIONAL RESULTS 

As mentioned earlier, the computation of the cohesive 
energies for various Fe, Mg, and also Si clusters in the 
bulk and at the Ga-terminated GaN(0001) surface allows 



CD 
it 

c 
o 

"-t— ■ 

-^ 
c 

CD 
O 

c 
o 
o 



+ 

CO 



CO 

'c 

Z3 
_Q 

a: 
</> 

LU 



(Ga,Fe)N 
(Ga,Fe)N:5Mg 



W H w X K 




4.3 4.5 4.7 4.9 

Field (kOe) 



0.11 
0.09 
0.07 
0.05 



• (Ga,Fe)N 

* (Ga,Fe)N:5Mg 



a 




50 100 150 200 250 300 350 
Cp 2 Fe flow rate (seem) 



FIG. 12: (Color online) ESR for (Ga,Fe)N:<5Mg (triangles) 
and reference (Ga,Fe)N (circles). Upper panel - selected Fe + 
transition and double integral for a digitally Mg-codoped (tri- 
angles, dashed line) and for a reference (Ga,Fe)N (circles, solid 
line) grown with 300 seem CP2FE flow rate. Lower panel - 
Fe 3+ concentration as determined by ESR (filled symbols) 
and SQUID (open symbols). 



us to interpret the effect of codoping by either Mg or Si on 
the incorporation of the Fe ions, especially in the frame 
of the different codoping modes. 



imately 1 nm. The first GaN bilayer is fixed in the ap- 
propriate bulk configuration. The unit cell is taken to 
be 2y / 3a x 3a, and the dangling bonds on the opposite 
nitrogen surface are saturated with 12 H atoms. The in- 
tegration over the Brillouin zone (BZ) is performed using 
a 4 x 4 x 4 Monkhorst-Pack k-point mesh. We have used 
the Methfessel-Paxton smearing for the integration over 
the BZ with a smearing width of 0.05 and 0.005 for the 
surface and bulk calculations, respectively. The atomic 
positions within the supercells have been fully optimized 
in all cases. The lattice constants are taken from our 
previous work^ and set to 0.4549 nm for zb-GaN, and 
a = 0.3211 nm and c — 0.5231 nm for wz-GaN. 

We have studied the formation of Fe„Mg m and Fe„Si m 
(n, m < 2) complexes in the zb-GaN matrix and also 
at the (0001) wz-GaN surface. The atoms in the stud- 
ied clusters occupy substitutional positions in the cation 
sublattice of GaN. The most stable configuration for the 
Fe-X (X = Mg or Si) cluster is the one where Mg or Si is 
the nearest neighbor (NN) of the Fe atom. If a second Mg 
or Si is at the NN distance to Fe, then the three atoms 
form an almost linear X-Fe-X cluster, where X = Mg or 
Si. A cluster with two Fe has been previously considered^ 
and the stable configuration has the two transition metal 
atoms as NNs, both in the bulk and at the surface. The 
Fe2X (X = Mg or Si) complex forms an isosceles triangle. 
Interestingly, the same configuration has been recently 
found to be the most stable structure for the isolated 
Fe 2 Si cluster ZL Finally, the Fe 2 X 2 (X = Mg or Si) cluster 
forms two isosceles triangles with a common Fe-Fe base. 
The two triangles are almost coplanar and the cluster 
has Cs symmetry. It should be pointed out that there is 
no significant difference in the shape of the structures in 
bulk compared to those at the surface. 



A. Methodology 



B. Results of ab initio studies 



The computations have been performed within the 
density functional theory-generalized gradient approxi- 
mation DFT-GGA+ U approach with the Perdcw-Burkc- 
Ernzerhof exchange and correlation functional 2 -* and ul- 
trasoft pseudopotentials^ to describe electron-ion inter- 
actions as implemented in the plane wave basis Quan- 
tum Espresso code<2& The energy cut-off is set to 30 and 
180 Ry for the kinetic energy and electronic density, re- 
spectively. The value of the U parameter is chosen to 
be 4.3 eV for the Fe 3d orbitals, a similar value being 
employed in our previous worki^ 

As previously^, we use the supercell approach to cal- 
culate the formation energy of Fe„Mg m and Fe„Si m 
(n, m < 2) clusters in the zinc-blende (zb) GaN and 
at the Ga terminated (0001) wurtzite (wz) GaN sur- 
face. The zb supercell consists of 64 (2a x 2a x 2a) 
atoms. A 4 x 4 x 4 Monkhorst-Pack k-point sampling 
mesh is employed. The modeling at the Ga surface is 
carried out with a supercell containing 4 GaN bilayers 
(48 Ga and 48 N atoms) and a vacuum region of approx- 



As previously reported^ two Fe ions tend to form pairs 
both in the bulk and at the Ga terminated surface as the 
heat of reaction (pairing energy E^) for those clusters 
is -0.158 and -0.125 eV, respectively. The formation of 
Fe-Mg and Fe-Si cation pairs in bulk GaN is energeti- 
cally favorable since E^ for those clusters is -0.178 and - 
1.349 eV, respectively. In contrast, at the Ga surface only 
Fe and Si tend to form pairs (E^ = -0.309 eV), whereas Fe 
and Mg avoid being at NNs positions (E^ = 0.187eV). 
It is worth mentioning that the E^ value for Fe-Si pairs 
in bulk is similar to the binding energy -1.44 eV for an 
isolated Fe-Si cluster. 27 Finally, neither Mg nor Si are 
likely to form monoatomic pairs in the bulk {E& =0.102 
and 0.847 eV, respectively) or at the surface (E^ — 0.130 
and 0.067 eV, respectively). 

The formation of Fe pairs may be, however, strongly 
affected by the presence of Mg or Si atoms. This is shown 
in Fig.1131 where we compare the heat of reactions for 
the formation of several consecutive Fe„X m (n,m < 2) 
clusters, where X denotes Mg and Si, respectively. 
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FIG. 13: (Color online) Heat of reaction (pairing energy Ed) 
for the formation (at negative Ed) of successive Fe n Mg m (up- 
per panel) and Fe n Si m (lower panel) nearest neighbor cation 
clusters (n, m < 2) in bulk zb-GaN (squares) and on the Ga- 
tcrminated (0001) wz-GaN surface (triangles). The solid lines 
are guide to the eye. The data points not connected by the 
solid lines (reversed triangles for surfaces, and diamonds for 
bulk) denote the heat of reaction of two nearest neighbor Fe 
cations in the presence of Mg or Si at the distance of 0.6 nm. 



Codoping of (Ga,Fe)N with Si or Mg not only affects 
the formation of the Fe pairs, but also influences the mag- 
netic coupling between the Fe ions, as shown in Fig. 1141 
where we plot the total energy difference between the 
ferromagnetic (FM) and antiferromagnetic (AF) config- 
urations for several Fe2X m (m < 2, and X denotes Si or 
Mg) clusters in the zb-GaN matrix and at the (0001) wz 
GaN surface. 

To complement our investigation, we have also car- 
ried out calculations for the formation of Fe n Mg (n < 
2) clusters in the presence of negative charge in the 
GaN matrix. The negative charge in the system may 
be viewed as an analog of the compensating defects 
{e.g., H impurities) that can be present during the 
codoping process of (Ga,Fe)N with Mg. The heat 
of reaction Fe + Mg~ — > FeMg~ is positive both in the 
bulk and at the surface (2?d = 0.012 and 0.196 eV for 
bulk and surface, respectively). The heat of the reac- 
tions Fe 2 + Mg" -> Fe 2 Mg" and Fe + FeMg~ -> Fe 2 Mg~, 
is 0.054 and -0.126 eV, respectively, in the bulk, and 0.514 
and 0.194eV, at the surface. This means that there is no 
qualitative difference between the results with and with- 



> 




UJ 



0.0 



-0.1 



Uf -0.2 



-0.3 




Fe, Fe Si Fe Mg Fe 2 Si 2 Fe Mg 



FIG. 14: (Color online) Energy difference between ferromag- 
netic and antiferromagnetic configurations for several Fe2X m 
(m < 2, and X denotes Mg or Si) clusters in zb-GaN and at 
the Ga-terminated (0001) wz-GaN surface. The solid lines 
are guides to the eye. 



out charge at the surface. However, in the bulk a neg- 
atively charged Mg is not likely to be attached to Fe or 
Fe 2 . For the hypothetical Fe 2 Mg~ the tendency to AF 
coupling between the Fe ions is preserved both in the 
bulk and at the surface. 

Finally, we have also studied the possibility of the for- 
mation of Fe pairs at the Ga surface in the presence of 
Si or Mg in the proximity of a Fe ion, in our case at 
a distance of 0.6 nm. The Fe ions - that have a nega- 
tive pairing energy at the clean surface - avoid pairing in 
the presence of Si or Mg, since E<± = 0.155 and 0.267 eV, 
respectively. The AF coupling is still dominant for the 
hypothetical Fe pair in the presence of Si or Mg at the 
surface. 



V. DISCUSSION AND CONCLUSIONS 

In this Section we discuss the outcome of the experi- 
mental studies summarized in Sec. Ill in the light of the 
theoretical results presented in Sec. IV. 

It is well known that except for isoelectronic Mn in 
II- VI compounds, where d-levels are far from the Fermi 
energy, in virtually all other transition metal (TM)-doped 
semiconductors a significant contribution of open d shells 
to the bonding leads to a strong attractive force between 
magnetic cationsi2i2£ This leads to segregation and/or 
precipitation of TM-rich compounds during the growth 
under thermal equilibrium conditions, the resulting sol- 
ubility limit being typically well below 0.1%. However, 
epitaxy at sufficiently low temperatures is known to al- 
low for the fabrication of dilute magnetic semiconduc- 
tors (DMSs) with random distribution of TM cations of 
concentrations largely surpassing the solubility limit, a 
canonical example being (Ga,Mn)As^2 The high tem- 
perature annealing of this material leads to the aggre- 
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gation of Mn cations resulting in the formation of Mn- 
rich embedded ferromagnetic nano crystals^ The ferro- 
magnetic/semiconductor nanocomposite formed in this 
way shows the typical properties of a superparamagnetic 
systems^ 

Extensive investigations over the recent years have 
made it increasingly clear that in most TM-doped 
semiconductors and oxides the aggregation of magnetic 
cations occurs already during the growth process.— The 
studies carried out to-date on TM-doped nitrides allow 
us to shed light onto a number of questions pertinent to 
these systems. 



trix, with the obvious effect of promoting the dilution of 
the system. 

It is worth noting that according to the results col- 
lected in Fig. [51 the presence of Fe favors the incorpora- 
tion of Mg into GaN. This may mean that the strength 
of the repulsive interaction within Mg-Mg pairs is low- 
ered by the presence of Fe. The ab initio results reported 
in Sec. IVB suggest that the pairing energies Mg-Mg and 
Mg-Fe are of the same order at the growth surface, 0.13 
and 0.2 eV, respectively. However, the presence of Fe 
leads to an attractive interaction within Mg-Mg and Mg- 
Fe pairs, that may promote the Mg incorporation, as 
observed. 



A. Does the aggregation of magnetic ions occur in 
the bulk or at the surface? 

A strong dependence of the efficiency of Fe aggregation 
on the growth rate led us to conclude that this process 
takes place at the growing surfaced This fact was fur- 
ther supported by the lack of evidences for aggregation 
of both Fe and Mn under post-growth annealing up to 
above growth temperature . 16 i 20 At the same time, we 
found that under similar growth conditions, the solubil- 
ity limit of Mn in GaN is at least one order of magnitude 
greater than the one of Fe^ We explained this surprising 
difference within first principles computations 9 implying 
that there is no attractive force between Mn cation pairs 
at the GaN surface, whereas the magnitude of the pairing 
energy of Fe cations E& = — 0.125 eV, reported in Fig. [13] 
can explain the formation of the Fe-rich nanocrystals, 
even at a growth temperature of ~1100K. 



B. Why does Mg, in contrast to Si, reduce the 
efficiency of Fe incorporation into the host? 

According to computational results displayed in 
Fig.1131 there is a gain in energy with the formation of 
Fe-Si cation pairs at the GaN surface. In contrast, there 
is a repulsive force between Fe and Mg under the same 
conditions. This means, in line with the experimental 
results collected in Fig. [5] at low Fe flow rate, that the 
presence of Si enforces the Fe incorporation whereas Mg 
has an opposite effect, as observed. 

At the same time, the destructive effect of homoge- 
neous Si and Mg codoping upon the Fe aggregation found 
experimentally in the previous 3 and present work - as ev- 
idenced in Fig. 0- is in accord with the theoretically es- 
tablished repulsive interaction between Fe cations at the 
GaN surface in the presence of either Mg or Si, as shown 
in Fig.[l3l As discussed in Sec. IVB, the effect is stable 
against a donor compensation {e.g., of Mg by H) and is 
independent of the distance of Mg or Si to the Fe-Fe pair. 
Once more it must be underlined, that in contrast to Si 
- shown to rather promote the incorporation of Fe into 
GaN - the homogeneous codoping with Mg diminishes 
the efficiency of the Fe incorporation into the GaN ma- 



C. How does digital doping promote the formation 
of nanocrystals? 

For the case of digitally codoped layers: (i) in con- 
trast to the case of homogeneous doping, the S- fashion 
growth-mode promotes the aggregation of Fe and the 
consequent phase-separation of the system; (ii) a high 
density of Fe-rich nanocrystals is observed already for 
0.1% of Fe under our growth conditions; (iii) the variety 
of the Fe-rich phases becomes more variegated with in- 
creasing Fe concentration; (iv) the nanocrystals are not 
distributed uniformly in the sample volume, but concen- 
trate 300 nm below the sample surface and at the in- 
terface GaN-buffer/(Ga,Fe)N:JMg-layer; and (v) due to 
the presence of the Fe-rich embedded nanocrystals, a sig- 
nificant FM contribution to the overall magnetization is 
observed up to above room temperature. At the same 
time, the presence of Fe reduces the mean concentration 
of Mg that can be incorporated in the digital mode. 

The origin of the above striking phenomena is not yet 
fully understood. It is possible that interdiffusion of Mg 
and Fe accompanying the digital growth mode makes 
that bulk rather than surface energetics applies to this 
case. According to Fig.[T31 the presence of Mg not only 
does not hamper the Fe aggregation but enhances it, par- 
ticularly in the case of larger clusters. At the same time a 
substantial repulsion between Fe-Mg pairs at the surface 
may lead to a reduction of a mean Mg concentration in 
(Ga,Fe)N, compared to digitally-Mg doped GaN. 

The overall effect of digital-doping and or/generally 
digital-growth requires at this point a dedicated study of 
growth-kinematics and -dynamics both in the presence 
and absence of intentional doping. 



D. Why are the embedded Fe-rich nanocrystals 
not distributed uniformly over the film volume? 

According to the present and previous 3 XRD and 
TEM studies, the average size of the embedded Fe-rich 
nanocrystals is independent of the Fe flow rate, sug- 
gesting that the precipitation occurs by a nucleation 
mechanism in which only nanocrystals with a critical 
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size can form. If this is the case, during the epitaxy, 
as long as the nanocrystals maintain an undercritical 
size, they develop and segregate with the growth front, 
a process that is interrupted once the critical size is 
reached. Accordingly, a high density of nanocrystals is 
expected to be found in narrow volumes perpendicular 
to the growth direction. Within this scenario, the 
distance between the nanocrystals-rich volumes would 
be inversely proportional to the rate of Fe incorporation 
at the surface. The previous and present TEM data 
reveal the existence of such 2D-like distribution of 
nanocrystals. An important consequence of the result- 
ing highly non-uniform distribution of the magnetic 
nanostructures over the film volume, is the onset of 
significant dipole-dipole interactions freezing the thermal 
fluctuations of the nanocrystals magnetic moments. 



VI. SUMMARY AND OUTLOOK 



The ab initio calculations, confirm the tendency of the 
Fe ions to aggregate at the growing sample surface of 
wz-GaN efficiently in the presence of Mg. Furthermore, 
they indicate that: (i) the presence of either Mg hinders, 
while the one of Si promotes the Fe incorporation, as ob- 
served; (ii) both (homogeneous doping with) Mg or Si 
quenches the aggregation of Fe. This is in contrast with 
the computational results for bulk, where an enhance- 
ment of the tendency to aggregation is expected in the 
case of Mg codoping. It is possible that this case ap- 
plies to the digital-like codoping, where interdiffusion is 
important. 

In general terms, our combined experimental and the- 
oretical studies demonstrate how particular codoping 
modes may serve to control the bottom-up fabrication of 
magnetic nanostructures in a semiconductor matrix. In 
view of our results, a considerable diversity of local atom 
arrangements obtained via particular codoping modes re- 
sults from substantial differences in the energetics at the 
growth surface and in the bulk, providing a worthy mean 
to design self-assembling processes at the nanoscale. 



The effect of Mg-doping on the structural arrangement 
and in the magnetic behavior of (Ga,Fe)N has been in- 
vestigated, with particular attention to the role of the 
doping method (homogeneous or digital). The onset of 
a competitive mechanism between Mg and Fe for the oc- 
cupation of substitutional sites becomes evident, as the 
presence of Mg reduces the amount of Fe in the homo- 
geneously Mg-doped layers compared to the correspond- 
ing reference samples (Ga,Fe)N and induces Fe aggrega- 
tion in the digitally fabricated ones. As a consequence of 
the reduced efficiency of Fe incorporation in the homo- 
geneously Mg-doped layers, the system becomes dilute 
and its magnetic response is purely paramagnetic. In 
contrast, upon JMg-doping, the tendency of the Fe ions 
to aggregate is enhanced, and the system becomes ferro- 
magnetic, due to the presence of FM £-Fe3N and galfenol 
Fei-zGa^ nanocrystals. 
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